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1 Introduction 

A large amount of effort has been denoted to 

capability - through the solution of elliptic 'forms near bound- 

present study, elliptic pdes used in gn gener . . diffusion-reaction problems, vortex 

aries that are similar to the equations used m nuclear “^^ition Mid convection) 

problems, electric space charge problems, s ea y 3 pde3 con tain appropriate inho- 

through long thin fins, etc. In the grid generation p oblems, these ptocwM ^ P.^ ^ ^ 

mogeneous terms that control the distribution o grx p referred to as Poisson equations 
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equations reduce to the long thin fin heat transfer equations with a finite heat transfer coefficient in 

the transverse direction (normal to the plane of paper) and a large heat transfer coefficient in the 
lateral direction. 

Most of the elliptic grid generation studies referred to above have been centered on developing: 
body conforming grids around bodies for external fluid flow simulations. The grids thus generated 
are smooth with at least first two derivatives continuous, appropriately stretched or clustered normal 
o any given coordinate direction and orthogonal over most of the grid domain. The inhomogeneous 
terms afford a grid control to satisfy clustering and orthogonality around specific surfaces (in three 
dimensions) and lines (in two dimensions). 

In external flows, these inhomogeneous terms, i.e., the source terms and the dependent variable 
proportional terms, are designed to vanish away from the body so the problem reduces to solving a 
Laplacian away from the body. 

In the present study, the inhomogeneous terms used are appropriate for an interior grid generation 
problem where all the boundaries enveloping the grid will affect the solution through these terms 
ihese terms are designed by interpreting their meaning physically through the principle of local 
conservation of thermal energy close to the grid boundaries. 

Two geometrical configurations, an annulus and a gear tooth, are considered to test the new 
oundary constraints derived in this study. The gear tooth geometry treated here corresponds to 
that ot a planar cross-section of a spiral-bevel pinion gear tooth typical of the OH-58 helicopter 
transmission pinion. This study is driven by the need to generate time-series vibration signatures 
rom the OH-58 helicopter transmission by finite difference simulation of the appropriate structural 
ynamic equations. The choice of elliptic pdes for grid generation is entailed by the need to generate 

ime series data as accurately as possible (see relative comparison with other representative grid 
generation methods in Ref. 7). ° 


2 Governing Equations 


The two-dimensional governing equations for an elliptic grid generation problem in an appropriately 
defined planar domain are 1 - 2 

txx + £yy = P (£, Tj) 

Vxx + Vyy = Q(£,T}) 

v are the generalized curvilinear coordinates, x and y are the Cartesian coordinates 
and the P(£ } r)) and are the inhomogeneous terms. 

The form of the inhomogeneous terms, P and Q, is, e.g., exponential 2 and is given by 


P(t,v) = 

- ai{T))sgn (f - £)exp(-6i|f - £|) 

(la) 

Q(i,v) = 

-Ci{Qsgn{T) - 7]i)exp(~di\T} - r]i\) 

(lb) 


where i refers to the grid boundary in question. 

For the sake of argument, without loss of generality, if we take the case where f > A and n > n- 
, then we have the inhomogeneous terms as " 


p (& = -cii(r))exp(-bi{€ - £)) 


(lc) 


<3(f , V) = -Ci(i)exp(-di(r] - 77*)) (id) 

At the boundaries, where £ = & and r, = m . Equations (lc) and (Id) respectively become 


l) = -a«(f?) 
and 

When 6;|£ — <f,| or d t jrj — rji j is small, the inhomogeneous terms take the form given by 

P&V) = -<*(t})(1 - bi(t - &)) 
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and 


Q(£,t l) = -Ct(0(! _ diiv ~ 

Therefore, the governing equation for, e.g., £, in the vicinity of the boundary £», becomes 

£*, + £ V * = — « i (^ 7)( 1 - &<(£ - 6 )) 


£ xx + iyy - o,i(rj)bi^ - -ai{rf) - ai(r])b t £i ( 2 ) 

If the term aA£, were absent, the resulting equation would turn out to be a Poisson equation. The 
equaticm 'given above arises, e.g., in the steady state heat conduction problems m long thm fins 
where £ is the temperature and where the heat transfer coefficient in the transverse thin dire t 
is moderate but is large in the lateral direction, and there is a balance amongst* he t heaU tonducted 
through the fin, heat carried away from or to it through convection m proportion to this moderate 
heat transfer coefficient and the heat sources/sinks distributed over the domain. Consider the case 
when £ > then, defining a new variable 


9 = £-£i 


Equation(2) becomes 


@xx *+■ @yy “ Cl{bi9 — CLi 


(3) 


The term can be interpreted as a heat source term. 

Equation(3) tells us that when £ > &, there is a balance between the heat converted from 
control volume in the interior to the boundary £ it heat conducted out of this control volume and 
heat lost from the control volume due to the heat sink, a* ( 77 V 

Similarly, if we consider the case when £i > £, then Equation (2) becomes 


e xx + Oyy - Oi(ri)bi0 = a* fa) 


(4) 


where 8 = & - £ 

The term a, can be interpreted as a heat sink term. 

Again Equation(4) tells us that when £ < there is a balance between the heat converted from 
the boundary £* to a control volume in the interior, heat conducted out of the control volume and 

the heat generated in the control volume due to the source, flifa)- . 

From SSons (3) and (4), it can be seen that for a given convective heat flux (given number 
of gridTiues), as the product, «A, decreases, the heat transfer coefficient decreases 

A v- u the temperature gradient at the boundary £ has increased so th 

?£££ ttZ ThL tii <J there is”a large gradient in f from the grid boundary i to 
the interior, thereby resulting in a highly entered l gnd j> the Laplace 

na,,y maintained there, since the 

“'tferrinTr^oMaTid (4), the Green's Theorem gives us respectively, 


1 j (-Oi + aibiO)da = d n 6ds 

(5) 

J J {a { + a,bi9)da = d n 6ds 

(6) 


where S is the surface area of a closed domain, C is the boundary enclosing this domain, n is me 
normal to the surface da is the elemental area and ds is an elemental arc. 

The integrands on the left hand side ±a t and o £ M represent the heat sink/source i term an^ 
the convection term respectively, and the integrand on the right han si e represen 
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solution process, which together with bi ensures the grid orthogonality and a given grid spacing at 
the ith grid boundary. 

In the design of the inhomogeneous terms, Pfarj) and Qfarj) in Eqs. 1(a) and 1(b), there is 
no restriction on the nature of the sink/source term, ia*. It can change sign which indicates the 
presence of sources and sinks, subject to the constraint given by Equations (5) and (6). Otherwise, 
improper combinations of sources and sinks will violate the extremum principle. If over the domain' 
there is a net rate of heat generation due to the source/sink combination, then there has to be a 
positive heat flux convected away and vice-versa. This requirement will automatically be satisfied 
by Equations (5) and (6). 

If there is a point heat source present in the domain, the isothermals (temperature or f or rj 
contour lines) will tend to cluster around it since the gradients in the vicinity of the source will be 
positive toward the source and high, depending upon the strength of the source, and conversely for 
a heat sink. Same argument applies to a line heat source and sink. By analogy, if the source term 
turns out to be positive over some parts of the domain, then the curvilinear coordinate lines will 
tend towards lines wdth higher coordinate values and vice-versa. 


3 Solution Procedure 

First, the boundary data are selected appropriate to the physics of the problem, so that the gradients 
in physical quantities can be resolved adequately. Since there is a symmetry plane and a rotational 
symmetry present in the gear problem, the half-tooth grid is reflected about this symmetry plane 
and then the tooth grid is rotated about a moving axis of periodicity, thus substantially reducing 
the computational effort in prescribing the initial distribution of the entire 19-tooth gear at time, 
t = 0. But, for the elliptic grid generation of the entire gear after this inital time, periodic boundary 
conditions are used in the circumferential direction. 

Then, by interchanging the independent and dependent variables, the governing equations to be 
solved in the computational space (£, 77) become 

~ + yXrm ~ — J 2 (P(£, rj)x£ Q(f, rfjx^) 

a VK ~ + iVm = + Q(€,v)y v ) 

These equations are solved in the computational space using a line SOR relaxation algorithm 
where each coordinate line in one curvilinear coordinate direction is solved semi-implicitly using 
the. Thomas algorithm for tri-diagonal systems. The inhomogeneous terms referred to above are 
designed and incorporated so that a desired grid behavior near the boundaries is achieved. 

The inhomogeneous problem is solved using a technique similar to that of Ref. 3 by over-relaxing 
the inhomogeneous terms during the iteration process. The inhomogeneous terms used in Ref. 3 
are well suited for external boundary value problems where they allow for clustering in only one 
curvilinear coordinate direction, normal to the body. But, in internal boundary value problems, 
inhomogeneous terms have to take account of the influence of the boundaries in both curvilinear 
coordinate directions. The inhomogeneous terms used here allow for clustering in both coordinate 
directions. 

The inhomogeneous terms, and are evaluated at the boundaries in terms of the 

left hand side at each line relaxation sweep. Then outward from each boundary, the inhomogeneous 
terms are attenuated through an exponential function in each direction, as discussed above. In f 
direction, outward from a given & boundary, this exponential term is of the form, -aiexp(-6i|£-&|), 
and, in rj direction, it is of the form, — aiexp(—bi\r} - 77* J). 

In what is stated above, a proof of concept study for an annulus and then a gear tooth is carried 
out to demonstrate the usefulness of the new boundary constraints. Then, the grid for the complete 
19-tooth gear is generated using periodic boundary conditions in the circumferential direction. This 
is the only way to solve the grid generation problem in a larger context of structural dynamical 
simulation of the gear, since the grid for the entire gear will be subjected to dynamic stresses 
nonuniformly. 

The boundary constraints given by Equations (5) and (6) are applied to a region close to the 
boundary. For example, a finite slender strip close to the boundary, 7? max , as shown in Fig. 1, is used 
as a control volume to evaluate the heat source term and the convective flux term over it with the net 
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, . hmiridarv the heat transfer coefficient, hi is large, 

heat flux calculated around it. At the n /max ! . fin is thin in t h e z direction, and the 

whereas the heat transfer coefficient, hi is moderate^ T^p^ ^ boundary can be imagined 

ambient enveloping the fin in the smaller, the temperature 

to be at a temperature, r ] - Vmax,*™ '- e - _ inrrpaaes for a mven heat flux (given number of 
gradient in the y direction near the boundary W* “he’ clustering there becomes denser. This 

grid lines, per the boundary prescription) an boundary becomes denser when 

is the physical basis for the argument that the clustering at a D i t - 

h La) parameter, fc, associated with the term m am a —uses. 



Figure 1: 


A schematic showing the finite slender strip over 


which the boundary 


constraint is applied 


4 Results 

■ * „ ...■.lir-uted bv generating two-dimensional interior grids for an 

The new boundary constraints for ° an annu l us , clustered at two opposing 

annulus and a gear tooth. Figure - t c f clustering at/near the t) boundaries, e.g., is 

boundaries in, e.g., the ( direction. « - wtee AS is the spacing at the r, boundary 

o„ s 

t i" «>— '■"« bKomes denser 

that boundary. mrameter convergence rate corresponding to 

Figs. 3a and 3b show the 3 ’ c L d 3d show the corresponding rate for the 

the clustering parameter value of 0.1, an 8 , number of iterations and the abscissa the 

upper boundary. The ordinate for t ese g^^ ^ stopped when the relative convergence criterion 

value of the decay parameter. Tne ca . within lOthe average initial spacing at that 

that the average final spacing at any given boundary is within lOthe g 

boundary. tn a c ta ble value rapidly. Corresponding to this 

Again, as is seen, the decay parameter com . 0 D boundaries are shown in Figs. 4a 

STL™ Z ZL Of grid points afong the { ditec, ion between 

the top and bottom r? boundaries. ,. . th c i us tering parameter value of 0-1 is 

— *> this 

grid. The corresponding decay function distribution is s own parame ter value of 0.05 is 

In Fig. 8, a grid for the gear tooth con^£S £ the convergence . Figure 

shown. Figs. 9a to 9d show the cor^ 0: | functio n distribu tion for the lower and upper 

10 shows the corresponding results tor tti. - -■ 

boundaries. 


o 
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Figure 2: Finite difference grid of annulus, clustering parameter = 0.1 



(c) 100 iterations (d) 1000 iterations 


Figure 3; Results for the annulus: (a), (b) Lower boundary decay parameter convergence, (c), (d) 
Upper boundary decay parameter convergence. 
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Figure 4: Results fo, the annulus: a) Lower boundary decay function distribution, b) Upper bound 
ary decay function distribution. 



Figure 5: Finite difference grid of pinion gear tooth, clustering parameter 0.1 
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(a) 100 iterations 





Figure 6. Results for the gear tooth, clustering parameter of 0.1: (a), (b) Lower boundary decay 
parameter convergence, (c), (d) Upper boundary decay parameter convergence. 
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Figure /. Results for the gear tooth, clustering parameter of 0.1: a) Lower boundary decay function 
distribution, b) Upper boundary decay function distribution. 
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Figure 8: Finite difference grid of pinion gear tooth, clustering parameter 0.05 



(a) 100 iterations (b) 700 iterations 



(c) 100 iterations ( d ) 700 lteratl0ns 


Figure 9: Results for the gear tooth, clustering parameter of 0.05: (a), (b) Lower boundary decay 
parameter convergence, (c), (d) Upper boundary decay parameter convergence. 
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Figure 10: Results for the gear tooth, clustering parameter of 0.05: a) Lower boundary decay 
function distribution, b) Lpper boundary decay function distribution. 


As c; f be observed by comparing with Figs. 5 and 6 corresponding to clustering parameter of 
‘ ’ t T 5 degree of clustering has increased at/near the boundaries, and therefore the decay parameter 

va ue has also correspondingly decreased when the clustering parameter has decreased from 0 1 to 
O.Od. 

Figure 1 1 shows a section of the complete 19-tooth gear grid which was generated using periodic 

boundary conditions in the £ direction (circumferential direction) with a clustering parameter of 
U.UD. 



Figure 11: Section of a finite difference grid of pinion gear, clustering parameter = 0.05 

The grid is seen to be smooth and orthogonal throughout in this case, as expected Figure 12 
shows the gear tooth grid for one-sided clustering at the Vmax boundary. The value of the clustering 
parameter at the lower boundary is 5.0, which means that the grid is coarsest at the lower boundary 
and it progressively gets refined toward the upper boundary. Dirichlet boundary conditions were 
applied at the & and £ max boundaries. The strict enforcement of orthogonality at these boundaries 
is not of concern here since the grid for the 19-tooth gear is generated with the periodic boundary 
conditions prescribed at these boundaries. Fig. 12 is another test case showing the efficacy of the 
new constraints derived here. 
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Figure 12: Finite difference grid of pinion gear tooth, clustering parameter _ 5.0 

Figures 13a and 13b show the corresponding lower boundary decay parameter convergence rate 
and Figs. 13c and 13d show the corresponding rate for the upper boundary. Here also, the decay 

Par Cmr2poTd V ing S to solution, the decay functions for the lower and upper boundarl “ 

are shown in Figs 14a and 14b. The difference here from the previous two-sided clustering is that 
the clustering at lower boundary, m , is sparse and that at the : upper boundary, 15 denSG 

desired. Accordingly, the decay parameter at the lower boundary is large and that at the upp 
boundary is small. The new boundary constraints in this case also yield a grid proper y c us r 

cit both the boundtirieb. m^dp with the 

4nother grid simulation with the same clustering parameter value of 5.0 was made with the 

complete 19-tooth gear. Periodic boundary conditions in the circumferential direction were used m 
this case. The corresponding grid is shown in Fig. 15. The slit .4.4' represents the radial line where 
periodic boundary conditions are enforced in the circumferential direction. 

4s a comparative example, a test case was run with an a priori prescribed value of the decay 

parameter of 0.15 corresponding to the case shown in Fig. 12, i.e., 

1 t-he resulting grid rapidly degrades as shown in big. lb(a) and the corresponding 

iTion in Jig. W). where after 50 iterations, the so, n, ion aborts 

due to negative 6 values of Jacobians. This is just an example of the problems encountered in the 
trial-and-error process in prescribing the decay parameter manually, for w ic now a so u ion 
been found through the new constraints derived in the present study. 

Finally, figure 17 shows a finite-difference grid model of pinion and driven gears in mesh. 
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(b) 1800 iterations 



Figure 13: Results for the gear tooth, lower boundary clustering parameter of 5.0: (a), (b) Lower 
boundary decay parameter convergence, (c), (d) Upper boundary decay parameter convergence. 




Figure 14: Results for the gear tooth, lower boundary clustering parameter of 5.0: a) Lower boundary 
decay function distribution, b) Upper boundary decay function distribution. 
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Figure 15: Finite difference grid of a section of pinion gear, clustering parameter - 5.0 




Figure 16: Comparative results for the gear tooth with a fixed decay parameter value of 0.1 
clustering arameter of 5.0: a) grid for the tooth b) convergence history. 
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Figure If : Finite difference grid model of pinion and driven gears in mesh 


5 Concluding Remarks 

The boundary constraints for elliptic grid generation problems developed in this study have been 
demonstrated for two internal geometrical configurations. Internal grids for an annulus, a gear tooth 
and a complete 19- tooth gear were generated using the new constraints. Smooth clustered grids have 
been generated using these constraints without hit-and-trial prescription of decay parameters and 
without any recourse to redistribution of grid points, which has been a common approach used in 
elliptic grid generation problems until now. With new constraints, elliptic grids can be generated 
in simulation time without any manual intervention thus making problems of structural dynamics 
and fluid dynamics over compliant boundaries straightforward. Thus, a fully automated elliptic grid 
generation methodology has been developed and validated. 
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